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Abstract. We study zeroes of Gaussian analytic functions in a strip in 
the complex plane, with translation-invariant distribution. We prove that 
the horizontal limiting measure of the zeroes exists almost surely, and that 
it is non-random if and only if the spectral measure is continuous (or de- 
generate). In this case, the limiting measure is computed in terms of the 
spectral measure. We compare the behavior with Gaussian analytic func- 
tion with symmetry around the real axis. These results extend a work by 
Norbert Wiener. 

Keywords: zeroes of Gaussian analytic functions, translation-invariance, 
ergodicity. 

1. Introduction 

Following Wiener, we study zeroes of Gaussian analytic functions in a strip 
in the complex plane, with translation-invariant distribution. Under certain 
assumptions on the spectral measure, Wiener proved that the zeroes obey the 
law of large numbers, and computed their horizontal density (limiting mea- 
sure). This result appears in his classical treatise with Paley [T71 chapter X]. 
Wiener's proof is quite intricate; this may be why it attracted little attention. 

In this work, we simplify Wiener's arguments and remove unnecessary as- 
sumptions on the spectral measure. We incorporate the result into a theorem 
that guarantees the existence of the horizontal limiting measure in question, 
and asserts it is not random if and only if the spectral measure is continuous 
or consists of a single atom. Then we prove a counterpart of this theorem for 
a natural class of Gaussian analytic functions which have a symmetry with 
respect to the real axis. 

For this purpose, we developed a general Edelman-Kostlan-type formula 
for computing the average zero-counting measure of zeroes of a symmetric 
Gaussian analytic function in some domain (see theorem [3] below). This result 
extends those of Shepp and Vanderbei [TH] , Prosen [TU] and Macdonald [13] ■ 
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1.1. Gaussian analytic functions. We will deal with two classes of random 
Gaussian analytic functions. 

Definition 1.1. Let D C C be a domain, and let {</> n }neN be analytic functions 
in D such that the series Yl n |0n(-z)| 2 converges uniformly on compact subsets 
of P. 

(1) Let a n be independent standard complex Gaussian random variables 
(a n ~ A/c(0, 1)). The random series J^An^n^) is called a Gaussian 
Analytic Function (GAF, for short). 

(2) Let b n be independent standard real Gaussian variables (b n ~ A/r(0, 1) ). 
If the domain D and the functions <p n are symmetric w.r.t. the real 
axis (the latter means that 4>(z) = 4>(z), z G D) then the random series 
Y^ n °n ( t ) n{.z) is called a symmetric Gaussian Analytic Function. 

Our assumptions on {4> n } ensure that the sums above a.s. converge to an 
analytic function in D j^l Chapter 2]. Throughout the paper we assume that 
there is no Zq G D such that 4> n {zo) = for all n G N (hence the function / 
has no deterministic zeroes). 

The covariance kernel of f(z) is defined by 
(1) K(z, w) = E(/( Z )/H) = E <f>n{z)JJw) ■ 

n 

The function K(z, w) is positive definite, analytic in z, anti-analytic in w, 
and obeys the law K(z,w) = K(w,z). It turns out that every such function 
K(z,w) of two variables z,w G D uniquely defines a GAF in D. 

If in addition K(x, y) is real whenever x, y G Dfll, then K(z, w) also 
uniquely defines a symmetric GAF with this kernel. We stress that a GAF 
and a symmetric GAF with the same kernel are different random processes. 

1.2. Stationarity. We assume our domain is the A-strip D = = {\lmz\ < 
A} with < A < oo. 

Definition 1.2. A GAF or symmetric GAF in a strip Da is called stationary, 
if it is distribution-invariant with respect to all horizontal shifts, i.e., for any 
t G R, any n6N, and any z±, . . . , z n G D, the random n-tuples 

(f( Zl ), . . . , f(z n )) and (f( Zl + t), . . . , f(z n + t)) 

have the same distribution. 

If f(z) is stationary in the A-strip, then for any x,y G R the covariance 
K(f(x)f(y)) depends on (x — y) only, so that K(x,y) = r(x — y) for some real- 
analytic function r : R — > C. From this we deduce K(z,w) = r(z — w) (both 
functions are analytic in z, anti-analytic in w, and coincide for z,w G R), and 
so r(t) has an analytic continuation to the 2 A-strip D 2 a- 



ZEROES OF TRANSLATION-INVARIANT GAFS 



3 



Since r(t) is continuous and positive-definite, it is the Fourier transform of 
a positive measure p (Bochner's Theorem): 



r(t)= / e 2mtA dp(X). 
Jr 

The measure p is called the spectral measure of the process f(z). 

Since r(t) has an analytical extension to the 2A-strip, p(A) has a finite 
exponential moment [121 Chapter 2]: 

/oo 
e 27r - 2Al|A| dp(A) < oo. 
-oo 

In fact, condition ^ is also sufficient for r(t) to have an analytic extension 
to the 2A-strip. Therefore, beginning with a finite positive measure p obeying 
pi), we can construct a kernel by 



(3) K{z,w) = / e 2wi{z -™ )x dp(\) 



which defines in its turn a stationary GAF in the A-strip. 

What measures could be spectral measures of a symmetric GAF? As we 
mentioned earlier, a kernel K(z, w) defines a symmetric GAF if and only if it 
is real for z, w G M; By relation ^ this is equivalent to p being symmetric 
with respect to the origin. 

Finally, we mention that a random GAF or symmetric GAF may be explic- 
itly constructed, as follows, from its spectral measure p. If {ip n (z)} n comprise 
an ortho normal basis in LVM.), then their Fourier transforms 



2nizX„ 



<j> n {z) = il> n (z) = / e ZmzX ip n (X)dp(X) 
Jr 

comprise a basis in the Hilbert space J-"{L 2 (]R)} (the Fourier image of L^iM.) 
with the scalar product transferred from L 2 (IR)). Then we have 

r(z-w) = E(f(z)J^)) = J2M* 



Therefore, when used in definition 1.1 the basis will give us a random 

function with the desired kernel. 

2. Results and Discussion 

2.1. Main Theorem. It will be convenient to introduce some notation: 

Notation 1. (zero-set, zero-counting measure) Let D C C be a region, and / 
a holomorphic function in D. Denote the zero-set of / (counted with multi- 
plicities) by Zf, and the zero-counting measure by n/, i.e., 

V0 E C™(D), [ <P(z)dn f (z) = V cf>{z) . 

Jd zez f 
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We use the abbreviation rif(B) = / drif(z) for the number of 

Jb 



zeroes m a 



Borel subset B C D. 



Notation 2. Let y € (—A, A). For a stationary GAF or symmetric- GAF in 
with kernel K(z,w), denote 

/oo 
e" 4 ^ A rfp(A) . 
-oo 

In the case of a GAF, define the function 

U W dy\tonP(3l)) d V [ f^e-^dpiX) J ■ 

In the case of a symmetric-GAF, define for y ^ the function 

(5) 



rfy ^tt^^) 2 - ^(0) 2 J dy 
and the positive number 



/ 



JZ Xe-^dp(X) 



Finally, a stationary GAF is degenerate if its spectral measure pf consists 
of exactly one atom. Similarly a stationary symmetric GAF is degenerate if 
Pf consists of two symmetric atoms (i.e., Pf — c(5 q + for some c, q > 0). 

The following theorem is our main result. Denote by mi the linear Lebesgue 
measure. 

Theorem 1. Let f be a stationary non-degenerate GAF or symmetric GAF, 
in the strip Da with < A < oo. Let v^t be a non-negative locally-finite 
random measure on (—A, A) defined by 



u f>T (Y) = ^n f ([0,T)xY), Fc(-A,A) 

Then: 

(i) Almost surely, the measures v^ T converge weakly to a measure v$ when 
T^oo. 

(ii) The measure Uf is not random (i.e. varz/j = 0) if and only if the spectral 
measure pf has no atoms. 
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(iii) If the measure Uf is not random, then: 

Vf = Lm\, if f is a GAF, 

Vf = S mi + R5o, if f is a symmetric- GAF, 

where 5q is the unit point measure at the origin. 

The measure Vf is referred to as "the horizontal limiting measure of the 
zeroes of /", or simply "the limiting measure". In the discussion and examples 
that follow, we assume the normalization ip(0) = j R dp(X) = 1. 

Remark 2.1. The limiting measure Uf might have atoms. Generally speaking, 
the weak convergence in the theorem guarantees that Vf^{[a,b)) converges 
to Vf([a,b)) for all a, b £ (A, A) with a possible exception of an at most 
countable set, which corresponds to atoms of the limiting measure ff. Yet, 
due to stationarity, in our case the limit exists on all intervals. We prove, as 
an example, the following result: 

Proposition 2.1. Almost surely, for any a, b £ (—A, A), we have: 

Jim v f , T {[a,b)) = u f ([a,b)). 

The proof is included in appendix |A} 

Remark 2.2. The part of the theorem pertaining to GAFs extends the afore- 
mentioned Wiener's theorem. In his work, Wiener assumed that the spectral 
measure p has the L 2 -density dp(X) = |0(A)| 2 <iA, that satisfies convergence 
conditions: 
For any \y\ < A, 



and 



(1 + x ) \(f)(x + iy)\ 2 dx < oo, 



'I + x 2 ) \ Q)'(x + iy)\ 2 dx < oo. 



As above, is the Fourier transform of <fi. Under these assumptions, Wiener 
proved that the limiting measure Vf exists and equals Lm\ where L is defined 
btl (§. 

Remark 2.3. (atomic spectral measure) 

Consider a spectral measure consisting of two atoms: 

The corresponding GAF is f(z) = (de" 2 " 92 + C 2 e 2 " 92 )/v / 2, where (1,(2 ~ 
Ac(0, 1), independently. The zeroes of such a function are 

c 2 



Zk 



Aug 



arg(^) + 2nk — i log 
Ci 



k £ N. 
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We see that all zeroes lie on the same (random) horizontal line, equally spaced 
upon it. The height of this horizontal line is a non-degenerate random variable, 
and so in this example vj is indeed random. 

For symmetric GAFs, the spectral measure above is degenerate (all zeroes 
of the corresponding function shall be real). We mention that it is possible to 
construct a random analytic function with continuous spectrum, for which a 
given asymptotic proportion of zeroes lie on the real line. For this, choose a 
continuous symmetric spectral measure, sufficiently close to the degenerated 
measure S q + 5- q . 

Remark 2.4. (behavior near the boundary and near the real line.) 

We observe that S(y) and L(y) have the same asymptotic behavior as y ap- 
proaches the boundary ±A. Therefore, zeroes of a GAF and of a symmetric 
GAF with the same kernel behave similarly near the boundary of the domain 
of definition. 

For a symmetric GAF, we observe a "contraction" of the zeroes to the real 
line: there are zeroes on the line itself, but they are scarce as we approach it 



from below or above (see figure 3.1 below). This is confirmed by a straightfor- 



ward computation, which shows that S(y) = 0(y), as y — > 0. 



2.2. Expected Zero-Counting Measures. In part(iii) of the theorem, the 
limit Vf(a, b) is actually the expectation En/([0, 1] x [a, b}). In order to calculate 
this quantity in the GAF case, we use the following classical formula, which 
appeared in Edelman and Kostlan's joint work on random polynomials [6]. 
Several proofs of this formula are known ([5], chapter 2). 

Theorem 2. (Edelman-Kostlan formula) For a Gaussian Analytic Function 
f with covariance kernel K(z,w), the expected zero-counting measure is given 
by 

(7) E(n f ) = ±AlogK{z,z). 

This should be understood as equality of measures in the following sense: 
for any compactly supported h G C°°(D), 

E J h(z)dn f (z) = ^- J h{z) A log K(z, z)dm 2 (z). 

Here and throughout this paper, m 2 denotes the planar Lebesgue measure. 

The proofs of this formula depend inherently on the fact that f(z) is a 
complex Gaussian random variable for all z, which fails for the symmetric 
GAF. To that end, we prove the following result, that extends previous results 
by Shepp and Vanderbei [19], Prosen [16] and Macdonald [13]. 
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Theorem 3. For a symmetric GAF f on some region with covariance kernel 
K(z,w), the expected zero-counting measure is given by 

(8) E(n/) = ^ A log (K(z, z) + y/K(z,z)* - \K( 

where the Laplacian is taken in the distribution sense. 

Notice that stationarity is not assumed in the last two theorems. Moreover, 
this formula combines information about real and complex zeroes. 

3. Examples 

3.1. Paley- Wiener Process (Sinc-kernel Process). Consider the spec- 
trum 

dp a {\) = — X[_ 0j0 ](A)dA, a > 0. 

Condition ^ holds for any A > 0, so the sample function / is entire. The 
kernel is: 

sm(2%a(z - w)) 

K{z, w) = — — — — = r(z - w) 

lixa^z — w) 



A base for construction of the GAF (in the sense of definition 1.1) is 



, , , sin(27ra2:) „ 
Ziraz — nir 

This example yields a surprising construction of a random series of simple 
fractions with known poles and stationary zeroes: Take for instance a — 1. 
Our function is 

sin(27T2;) 



/(*) = £■ 



2tcz — rnr 

where {a n } are independent Gaussian random variables. Almost surely, Zf D 
|Z = 0, so we may divide by sin(27r^)/7r and get the random series 

g( Z ) = \^ . 

The poles of g are known (and lie on a one-dimensional lattice), but its zeroes 
are a random set invariant to all horizontal shifts! 

Using Theorem [T] for p a , we get that the zero-counting measure has the 
following density of zeroes: 



- (£a) =4na2 aty{ C ° thy -l) 



Similarly, the symmetric GAF with the same spectral measure has the con- 
tinuous density of zeroes 



4naJ dy \ x / s [ n h 2 y - y* 



I a/ sinh 2 
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plus an atom at y = 0, of size 




-15 -10 -5 5 10 15 -10 -5 5 10 -0.2 -0.1 0.0 0.1 02 



(a) Paley- Wiener (b) Fock-Bargmann (c) exponential spectrum 

Figure 1. Horizontal density of zeroes for GAF and symmetric 
GAF models with the same kernel. In each model, the lower 
graph represents the continuous component of the mean zero 
counting-measure for the symmetric GAF (the atomic part is 
an atom at y = 0, which is not graphed). The upper graph 
represents the continuous (and only) part of this measure for 
the appropriate GAF. 



3.2. Fock-Bargmann Space (Gaussian Spectrum). Set 

dp a (\) = -^=e- x2/a2 d\ , a > 0. 
Once again, / is entire. The Fourier transform of the measure is 

therefore the covariance kernel is: 

K(z,w) = e -^^)\ 



This space has an orthonormal basis of the form 




2 , where b = v/2- 



and c = 

In this model, the density of zeroes is constant: 

L - (tt-) = 2 ™ 2 
\27ra/ 

This is the only model with Lebesgue measure as expected counting measure 
of zeroes. For more information about this model and why the distribution of 
zeroes determines the GAF, see [22], [IS], P3] or [5J Chapters 2.3, 2.5]. 
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However, for the real coefficients case the continuous part of the limiting 
measure has density 



S a 



y 

2na 



2na 2 



d 



dy \ Ve 2 y 2 - l 



and the atom at y — is of size 

Both continuous densities are graphed in 1(b) , for the parameter a 



4tt ' 



3.3. Exponential Spectrum. Consider a symmetric measure with exponen- 
tial decay: 

dp(X) = sech(7rA)<iA 



-dX 



cosh(7rA) 

then r(z) = sech(7rz) as well. This model is valid in the strip —j< Im(z) < \. 

m 



Here 



cos 2 (27ry) 

For the symmetric GAF in this model, we have 

?r| sin(27n/)| 
cos 2 (27iy) 

We see that the zeroes concentrate near the boundaries of the region of con- 



vergence (figure 1(c)). 



4. Proof of Theorem [3]- Zero-Counting Measure for a 

symmetric gaf 

In this section we prove Theorem [3] Similar formulas were proved in specific 
cases. Our proof follows Macdonald [IB] , who has considered random polyno- 
mials (also in mult i- dimensional case). A novelty is in extension of his result 
to arbitrary symmetric GAFs. 

Recall that for any analytic function / (not necessarily random) in a domain 
D we have 

1 



ri 



f 



2tt 



Alog|/|. 



This is understood in the distribution sense. 

Using this for our random /, we would like to take expectation of both sides, 
to get: 



E 



(9) 



h(z)dnf(z) 



x 



E 



1 



1 



Ah(z) \og\f(z)\dm 2 (z) 
Ah(z)E[log\f(z)\]dm 2 (z) 



where m denotes the Lebesgue measure in C. The last equality is justified by 
Fubini's Theorem, as we will show in the end of this section. Thus we can 
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conclude that (in the weak sense): 

(10) E(n / ) = i-AElog|/|. 

Let us return to our setup: / is a random function generated by a basis <pk{z) 
of holomorphic functions, each real on R, and such that the sum E& |(/>j,(z)| 2 
converges locally-uniformly. Denote <f>k(z) = Uk(z) + ivk(z) where Uk,Vk are 
real functions. Our random function is decomposed thus: 

f( z ) = ^2 b kM z ) = ^2 b kU k (z) + % ^ hvk(z) = u{z) + iv(z), 

where b^ ~ A/r(0, 1) are real Gaussian standard variables. (u(z),v(z)) have a 
joint Gaussian distribution, with mean (0,0) and covariance matrix 

E u l E u kVk 
E u kVk E v k 
Lemma 4.1. The matrix E from above has 2 eigenvalues: 

K(z,z)±\K(z,z)\ 

A 2,l — ~ 



where K(z,w) = X>fcW0fcM = E 4>k{z)4)k{w) . 
Proof. For any complex number = u + iv , we have: 



M 2 = I(|0| 2 + Re(0 2 )), v 2 = I(|0| 2 -Re(0 2 )), m; = ^lm(0 2 ). 
Applying this, we can rewrite E as 

y _ n(EW 2 + ReE0 2 fc ) |lmE0 2 fe A, 

V l Im E^ |(EI0 fc | 2 -ReE0l) 7' 

and then calculate it's determinant and trace: 

AiA 2 = det E = \ ((£ |0 fc | 2 ) 2 - (Re$> 2 )) 2 - (lm^(0 2 )) 2 ) 

(11) =l(^,z) 2 - \K(z,z)\) 

Ai + A 2 = trace E = ^ |0 fc | 2 = if (2, 2) 
The lemma follows. □ 
Using the law of bi-normal distribution, we get: 

(12) E[log|/(z)|] = j / f \og{^f^^)e-^ y) ^ 1{x ' y)T dxdy = 

27rvdet E J Jr2 

^=\l log^x 2 + y 2 )e-^ lx2+x ^ y2) dxdy 

2ir v det E J Jr2 

Applying to the last integral the change of variables x = u^/Xl, y = w-yA^, 
with Jacobian ^/A~iA~2 = Vdet E, we have 

E[log|/(z)|] = -^ ff MVAxu 2 + \ 2 w 2 )e^ u2+w2 ^ 2 dudw 

^ J JR2 
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Now, changing to polar coordinates u = rcos9, w = rsin# we get: 
Epog|/(z)| 



11 



- / log( cos 2 9 + A 2 sin 2 0) dQ + C, 
!7r Jo 



2vr 

where C is a constant which does not depend on the point z (i.e., is independent 
of Ai and A2). In the following, we write C for any such constant (which may 
be different each time we use this symbol). These constants will vanish when 
we apply Laplacian (recall (Jio|)). 

So, the integral we should compute is: 



2 k 



log 



Ai cos 9 + % a/ A2 sin 6 

2tt 



d6_ 
2^ 



rzir 

\og(>/>^+ a/A^) + / log 
Jo 



e 2ie + v Ai - vA 2 



Ai + VA 2 



dQ 
2^ 



C. 



The remaining integral is computed easily by Jensen's formula for the func- 
tion z 2 + c, where c = v ^P-7 v aP^ < 1- Indeed, it has two zeroes in the unit circle, 



denoted a% and a 2 , and so: 



2tt 



log \e 2ie + c| 



dd 
2^ 



log 



— log |ai| — log |a 2 | = log \c\ — 2 log y\c\ = 0. 



Recalling (10) and using the relations (11), we arrive at 

E ( dn f) = ^ I lo S ( A i + A 2 + V4A1A2) = Alog z) + v 7 ^^,^) 2 -!^ 

4.1. justification of (|9]). We must show that the following integral converges: 
1 



2tt 



A /i^)|-E| log |/(*)| \dm 2 (z). 



x 



It is enough to prove that E| log \f(z) \ | is bounded on a compact subset S 
of the plane. f(z) is a 2-dimensional real Gaussian variable with parameters 
noted above, so we get 

E|log|/(z)[[= L_ ff \\og(^TV 2 )\e-^ lx2+X * ly2) dxdy 
27rvdet L J Jm.2 

As before, Ai,A 2 are the eigenvalues of S, dependent on z. By another 
change of variables (x = wa/Ai, y = w\0^) we get: 

E| log |/(^)|| = — // |log(AiM 2 + X 2 w 2 )\e- {u2+w2)/2 dudw 
4vr J J R2 

Fix z, and assume Ai < A 2 . Let us split the integral into two domains: 
Q + = {(u,w) G R 2 : log(Ai« 2 + \ 2 w 2 ) > 0} and = {(u,w) E M 2 : 
log(A lM 2 + X 2 w 2 ) < 0}. 

Then, on fi + we estimate < log(AiW 2 + A 2 u> 2 ) < log(A 2 ) + log(w 2 + w 2 ). 
Fro m here clearly the integral on fl + is bounded by C + G\ log A 2 . By lemma 
4.1, A 2 = ^(K(z, z) + \K(z, z) |) is a continuous function of z, and therefore is 



bounded on our compact set S. 
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For (u, w) G ft- notice that > log(AiW 2 + X2W 2 ) > log(AiW 2 ), therefore: 




log(Aiu 2 + X 2 w 2 )\e- {u2+w2)/2 dudw < 



// (| log(Ai) + \og(u 2 )\)e-( u2+w2 V 2 dudw <C + CillogAil 

Denote m = min{Ai(z) : z G S}. If m = 0, this leads to K(zq,Zo) = for 
some zq G K, but this means zo is a deterministic zero. Therefore m > and 
I log Ai| is bounded from above. 

5. Proof of Theorem Q] - Horizontal Limiting Measure 

5.1. Preliminaries. We present the probability space of our interest, equipped 
with a measure-preserving transformation. We explain the notion of ergodicity 
in this setup. 

The probability space Q is a countable product of copies of C, with P being 
the product of complex Gaussian measures (one on each copy). These copies 
represent the random coefficients in the construction of /: each u = {a n } n G Q 
corresponds to a function f u (z) = ^a n 0„(2i). Tf is the Borel a-algebra 
generated by the basic sets {u G Q : f u (z) G B(w,r)}, where z G D,r > 0. 
Here B(w,r) = {p G C : \p — w\ < r}. The group of automorphisms 5^ shall 
be defined via the correspondence to -H- f^: 

fs t M) =fu(z + t). 

The map St is measure-preserving, since we assumed that / is stationary. 
Thus, we will say the random process f(z) is ergodic, if any measurable set 
A G Tf which is invariant to all translations {St A — A, Vt G K) is in fact 
trivial (FA G {0,1}). 

In a similar way, one can define when is the zero-set Zf ergodic (it is itself a 
random point-process in the plane). The space Q, the measure P on it and the 
automorphisms {S t } are just as before. Now, the a- algebra J^z f is generated 
by the basic sets {co G Q : Zf w n B(z, r) ^ 0} with z G D, r > 0, B(z, r) C D. 

Corollary 5.1. Ergodicity of f implies ergodicity of Zf. 

Proof. It is enough to prove Tz f C Tf. Let A be a countable dense set in C. 
Basic sets of Tz, can be written as 




which is indeed in J-}. □ 

We will use the following classical result: 

Theorem 4. (Fomin, Grenander, Maruyama) A stationary GAF (symmetric 
or not) is ergodic w.r.t. horizontal shifts {f(z) —> f(z + t)}} t£ R if and only if 
its spectral measure p has no atoms. 
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This theorem was originally proved for real processes over M. (see for instance 
Grenander [7]), but small modifications extend it to a strip in the complex 
plane for both types of functions (GAFs and symmetric GAFs). 

5.2. Existence of the horizontal limiting measure (statement (i) ). 

As above, for T > 1 let be the random locally-finite measure on (—A, A) 
defined by: 

(13) ^ = >„(y). n ^ ,r) , FC(-A,A). 

In this section we show that a.s. the measures converge weakly as T tends 
to infinity. First, we assume that T tends to infinity along positive integers. 

In this case, we use the subscript iV instead of T. By a known theorem in dis- 
tribution theory (see for instance, [8j section 2.1]), a sequence of measures 
converges weakly to some measure if and only if the sequence of real numbers 
VN{h) is convergent for every h G Cq°(— A, A). It suffices to check whether 
VN{h) is convergent for all h G M, where M C C£°(— A, A) is a dense set of 
test-functions, and we may choose M to be countable. Given test function 
h G M, denote by Ah the event that vn(K) is a convergent sequence of num- 
bers. To prove our claim it suffices to show P(^U) = 1 for every h G M. Note 
that is N (h) = ±(X 1 + X 2 -\ h X N ) where 



(14) X k = X k (h) = j l [k<k+1) (x)h(y) dn f (x } y) 

is a stationary sequence of random variables. 

The random variables X k are integrable. This follows at once from an 
Offord-type large deviations estimate (5], theorem 3.2.1]: 

Theorem 5 (Offord-type estimate). Let f be a Gaussian analytic function on 
a domain D. Then for any compact set K C D, the number rif(K) of zeroes of 
f on K has exponential tail: there exist positive constants C and c depending 
on the covariance function of f and on K such that, for each A > 1, 

F{n f (K) > A} < Ce~ cX . 

Therefore we can apply Birkhoff theorem [U chapter 7]. It yields that the 
limit jr(Xi + X 2 + ■ ■ ■ + X N ) almost surely exists, and so f(A h ) = 1. This 
completes the proof of the weak convergence of the sequence vjq. 

Now, we consider the general case in statement (i). Let T > 1, and let 
iV = [T] be the integer part of T. Then 

N If 

v T {h) = Y u N(h) + - / i[ NtT )(x)h(y) dn f (x,y) . 

v .. ' 



=:R T (h) 

We will show that a.s. the second term on the right-hand side converges to 
zero for all bounded compactly supported test functions h. It suffices to prove 
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this for all bounded test functions supported by an interval [— Ai, Ai] with an 
arbitrary Aj < A. We have 

\R T (h)\ < H^n/Qtf.tf +1) x [— Ai, Ai]) . 

Employing the Offord-type estimate with K = [0,1] x [— Ai,Ai] and using 
translation-invariance of the zero distribution of /, we see that for each e > 0, 

F{n f ([N,N+l)x[-A 1 ,A 1 )) > eT) = P{n/([0, l]x[-A 1} A 1 ]) >eT} < Ce~ ceN 

Hence, for each MgN, 



P{limsup| J R T (/i)| > e||/i||oo} < Y] Ce~ C£N = C(l - eT^e 
and we conclude that a.s. 



ceM 



lim R T (h) = 

T— s>oo 

for all smooth compactly supported test functions h. This completes the proof 
of statement (i) in theorem [T] □ 

5.3. Non-random limiting measure (statement (iia),(iii) ). Here we will 
prove that if the spectral measure pj has no atoms, then the horizontal mean 
zero-counting-measure i/f is not random, which is a half of statement (ii) (the 
case when pf is a single atom for GAFs or two symmetric atoms for symmetric 
GAFs was already discussed in remark 2.3). We then compute the limit Vf, 
which is statement (iii). 

Assume the spectral measure p is continuous. From the Fomin-Maruyama- 
Grenander Theorem we get that / is ergodic, and so is Zf. Using the notation 
introduced in the proof of statement (i), we get that for any smooth test 
function h, the stationary sequence of random variables Xk(h) introduced 
in (14) is ergodic. In this case the Birkhoff ergodic theorem asserts that, a.s., 

1 N ^ 

lim u N (h) = lim - V X k (h) = EX (h) . 

N^-oa N^oo iV ^ — ' 

k=0 

Therefore, the horizontal mean zero-counting-measure is non-random, and 
equals 



EX (h)=E j i [0A) (x)h(y)dn f (x,y) = j l [0A) (x)h(y) Edn f (x,y) , 

where Edrif(x, y) is the mean zero-counting measure. For a GAF, we compute 
it directly by Edelman-Kostlan formula ([7]); while for a symmetric GAF we use 

/■oo 

3). As before, denote ip(y) = / e" 47r ^rfp(A). 

J — oo 



the formula rt8|) (in Theorem 
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Note that 

K(z, z) = J e 2 ^ 2yt dp(t) = i[)(y) 

K{z,z) = [ e^-^dpit) = [ dp{t) = V(0), 



where z = x + iy. Putting this into (|8j), we get the first intensity of zeroes: 

En, = ^^rlog U{y) + ^(y) 2 -^(0) 2 ) = j-j- . ^ 
1 Andy 2 *V yyi VV ^ y; Air dy ^(y) 2 - ^(0) 2 

For any y ^ 0, this is a derivative in the functional sense, which equals S(y). 
At y = 0, the function is not defined; but the limits 

Atx ^Jip(y) 2 - ifj(0) 2 y^°- An ^(y) 2 - ip(0) 2 

exist. This follows from f M == being an odd function, increasing in 

V^(y) 2 -^(o) 2 

y G (0, A). So, in order to compute En, we take the required derivative in the 
distribution sense, which yields the continuous point-wise derivative S(y) (for 
y 7^ 0) plus an atom of size 2 A at 0. 

In order to compute A let us write this limit again, and apply L'Hopital's 
rule: 

*A = Urn iM = . im ^^WEgg = ( If . . 1 



A 2 rfp(A) 

where £2 = — ?H r~^~ ls ^ ne ra tio between the second and the zero moments 

J-00 rf P( A ) 

of the spectral measure. We conclude that A = \f&n and therefore the atom 
has twice this size. 

5.4. Random limiting measure (statement (iib)). Our goal is to prove 
the second half of (ii). We present the proof for symmetric GAFs, since it is 
slightly more involved. We assume that the spectral measure has the form 

Pf = c $q + c $-q + £*> 

where p is a non-trivial measure. Our goal is to show that the horizontal 
mean zero-counting measure of some segment Uf(a, b) is a non-constant random 
variable. We may assume that c = 1 and q = 1 (if q = the analysis is easier). 

Since L 2 {R) is the direct sum of L 2 ^^ (R) and L 2 (IR), a union of any 
orthonormal bases in these subspaces is an ortho normal basis in L 2 (IR). By 



the remark in end of section 1.2, applying Fourier Transform on this union 
we shall get a basis <j) n (z) from which a GAF with spectral measure p can be 
constructed. This gives the representation 

f(z) = g(z) + acos(27T2;) + (3 sin(27rz), 
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where g(z) is a symmetric GAF with spectral measure \i and a, {3 ~ A/r(0, 1) 
are real Gaussians, independent of each other and of g. We write for short 
r](z) = r] a ^(z) = acos(27T2) + (3sm(2nz). 

Fix a,b G (—A, A). Denote the number of zeroes of / in [0, T] x [a, b) by 
N T {g,a,P) = #{ze[Q,T\x[a,b):g(z) = -Ti 0h f i {z)}. 

Assume to the contrary that there is some constant C (depending on a and 
b) such that 

o 3 r N T (g,a,fi) 
(15) a.s. ma, p, 3 hm — = C. 

Here we denote P 9 and E s the probability and expectation (resp.) condi- 
tioned on a, (3. We claim that: 



(16) 



E„ Urn = lim WW) 



T 



This exchange is justified by the dominated convergence principle, as seen 
by the following Offord-type estimate: 

Proposition 5.1. Let g be a symmetric stationary GAF on a horizontal strip, 
a and (5 are fixed complex numbers. There exist positive constants C and c 
such that: 

supr J!Mi^l >s]< ce- 

T>1 V 

This fact is proved in section [B] below. 



Next we claim that the RHS of (16) is just K g Ni(g, a, /3). To see this, 
notice that for integer T, NT(g,<y,/3) is the sum of T identically distributed 
random variables, all distributed like Ni(g, a, /3). This follows immediately 
from stationarity of g and from 1-periodicity of T] a p(z). Therefore for integer 
T, 

^E g N T (g,a,P) = E B N 1 {g,a,P). 

Denote M — [T\. Since E g N [M ,T] ■= G [M, T] x [a, b) : g{z) = 

—Va,i3(z)} < E 9 A^i < oo, it follows that for non-integer T, 



17) ton MH^i) = lim (Mm . E 



(18) 



Combining (15), (16) and (17) we have: 

a.s. in a, j3, ~E g Ni(g,a, f3) = C. 



We divide the rest of our argument into three claims. 
Claim 5.1. K g Ni(g, a, (3) is continuous in (a, (3) G M 2 . 
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Claim 5.2. For any compact set K C D, let N(g,a, (3; K) be the number of 
solutions to g(z) = —rja^z) with z G K . Then 

E g N{g, a, (3; K) n cos(27Tz) (K) as \a\ oo. 

Here n cos ( 27r2 ) is the zero-counting measure of cos(27r z) . 



Relying on the last claim and (18), we get that 
(19) E g N 1 (g,a,P) = 25 ([a,b)) 



for almost all a,/3. Since EgN^g, a, (5) is continuous in a, f3, (19) is true for 
all a, (3, and in particular for (a, (3) = (0,0). The following claim asserts this 
happens only for one family of symmetric GAFs: 

Claim 5.3. If for -A < a < < b < A, 

E g N 1 (g,0,0)=En g ([0,l) x [a, b) = 25 Q ([a, b)), 

then the spectral measure of g is |(<5i + 5_i) (up to a constant multiplier) . 

From this last claim it follows that the spectral measure of / consists only 
of symmetric atoms at ±1, which contradicts our assumption. 

It remains now to prove the claims. In the course of their proof we justify 
exchange of limits and expectations by the following 

Proposition 5.2 (Offord-type estimate for sine-like levels). Let g be a sym- 
metric Gaussian analytic function on a domain D, and let a and (3 be fixed 
complex numbers. Then for any compact K C D, the number N(g,a, f3; K) 
of solutions to g(z) = —T} a> p{z) with z G K has exponential tail: There exist 
positive constants C and c such that 

F(N(g,a,P;K) > s) < Ce' 08 . 



This proposition's proof is much similar to that of proposition 5.1, and we 
omit it. 



Proof of claim 5.1 Fix —A < «o < A> < A. It is clear that almost surely in g, 
N\(g, a, Pq) approaches Ni(g, a , /3 ) as a approaches a (the event of having 
a solution to g(z) = —r] a0j/3o on the boundary of [0, 1] x [a, b] is negligible). 



By proposition 5.2, we may pass to the limit: 
lim E j; iV 1 (0,a,A>) = lim / ^(N^g, a, (3 ) > s)ds 



lim P^iVxfoar, A>) > s)ds = [ f g {N x {g, a , Po) > s)ds = E g (N 1 (g, Co)). 



□ 



Proof of claim Fix p and g. For any q^O, the zeroes of 

, . , q(z) + sin(27rz) 

hJz) = — + cos(2ttz) 

a 



18 



NAOMI FELDHEIM 



and of f(z) = g{z) + Va,p{ z ) are identical. Now notice that h a (z) converges 
locally uniformly to cos(27T2;) as a — > oo (i.e., uniformly on any compact 
set). By Hurwitz's Theorem, this implies that the zero counting measures also 
converge locally uniformly, in the sense that for any compact K C D, 

lim n ha (K) = n cos{2nz )(K). 



By the bound in proposition 5.2 this almost sure convergence in g yields 
moment convergence: 

E g n ha {K) -> n cos{27rz )(K), as a -> oo. 

□ 




R=2a / X 2 dp(X) 



Proof of claim 5.3 Suppose the spectral measure is normalized, so that tMO) = 
J R dp(\) = 1 (else, multiply it by a constant). The premise and theorem [llgive 
two conditions on ip{y) = K(iy,iy): 

for some constant c E R. Solving the left-side ordinary differential equation, 
and using -0(0) = 1, we get ip(y) = cosh(q/). Since ip is a Laplace transform of 
p, we get p = \{5 c /2tt + 5- c /2tt)- But the right-side equation would be satisfied 
only if c = 2n. □ 

Appendix A. Convergence on all intervals (proof of 

proposition 



2.1) 



In this section we prove proposition 2.1| We use the notations developed 
in section 5.1 For any point in the probability space u> G Q, let be the 



sequence of measures introduced in (13), for integer T = N (the non-integer 
case follows just as in the proof of part (i) of Theorem [TJ and will not be 
discussed). 
Define the set 

C = {u G fl : Wn)n converges weakly} 

Notice that by part (i) of TheoremjlJ P(C) = 1. For convenience of notation we 
denote the limit Vf = . From general measure theory, one can deduce that 
almost surely, u^Qa, b)) converges to z^([a, b)) for all a, b out of a countable 
exceptional set. This exceptional set is the set of atoms of vj (which might be 
random). We thus turn to define 

A = {u> G Q : lim z^ia} = z/^{a}, for each atom a of z/ w } C C. 

Claim A.l. A is measurable with respect to J 7 /. 

The proof of this claim will be presented in the end of this section. Our 
next goal would be to prove: 
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Claim A.2. F(A) = 1. 

Our main tool will be the Ergodic Decomposition Theorem (proved, for 
instance, in [TJ chapter 2.2.8]): 

Theorem 6 (Ergodic Decomposition). Let (Q,J-,¥) be a standard Borel- 
space, equipped with a measure preserving transformation S : Q — > Q. Then 
the set E S (Q) of ergodic probability measures on Q is not empty, and there 
exists a map (5 : Q — > E S (Q) such that for any measurable set A G J 7 the 
following holds: 

(1) the map < / A is measurable. 

(2) P(A) = / n /3 w (A)rfP(a;). 



Proof of claim A. 2 . The stationary system (f2, J 7 ^, P, 5) defined in section 5.1 
and the set A defined above meet the requirements of the Ergodic Decompo- 
sition Theorem. Therefore, in order to prove our claim it is enough to show 
that 

V 7 G E s (Q),-f(A) = 1. 

Fix an S'-ergodic measure 7. Since A is an ^-invariant set, we get 7(A) G 
{0, 1}. Moreover, the event has an atom in the interval /} is also invari- 
ant, for any interval / C (—A, A). Therefore, 7-a.s. the limiting measure 
has atoms at some known points (a n )„ e pj C (—A, A). 

For a certain atom a = a n , define the stationary sequence: 

X k (a) = n f ([k, k + 1) x {a}), 

and notice that 

v N {a} = j^^2 x k(a). 

As 7 is ergodic, we have by Birkohff's ergodic Theorem: 

7-a.s. v%{a} converges to E 7 n/([0, 1) x {a}) = ^{a}, as N — > 00 
Since there are at most countably many atoms, we get 'j(A) = 1. □ 



We now know that P-a.s., the sequence is weak convergent and converges 
on any atom of the limiting measure (to the desired limit). A general claim 
from measure theory will assure us that in this case, converge on any 
interval: 

Claim A. 3. Suppose {vn)n is a weak-converging sequence of measures on 
some interval I, and let v be the limiting measure. If lim = v{a\ for 

every atom a of v, then lim Vn{J) — ^(J) for every interval J G I. 
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Proof. We demonstrate the case J = [a,b), where v has no atom at b (other 
cases are similar). 

Given e > 0, one can construct piecewise linear functions + , 0~ G C{I) 
such that: 

(20) Wx, <fr{x) < 1 (0)6) < l [a>b) (x) < <f> + {x), 

and additionally 

< u((j>+) - {v(4T) + v{a}) < e. 

(For instance, for large enough parameter n, take (f) + supported on [a— -, b], 
equals 1 on [a, b — -] and linear otherwise; <p~ supported on [a, b], equals 1 on 
[a + -, b — -], and linear otherwise). 

By applying the measure to relation (20), we get: 

v N {4T) + v N {a} < u N ([a,b)) < vn{<P + ) 

But, from our assumptions, for large enough N we have 

u{(j>~) + v{a} -e<is N ([a, b)) < z/(0+) + e 

As the difference between those bounds does not exceed 3e, we see the limit 
lmi/v->-oo ^jv([0) b)) exists. Since ^(0 + ) is as close as we want to v([a,b)), we 
are done. 

□ 

It remains only to prove the measurability of A. 



Proof of claim \AJ\ We first investigate some underlying objects. Denote by 
P = P{— A, A) the space of all locally finite measures on (—A, A) induced with 
the Levy - Prokhorov metric (for which convergence in metric is equivalent to 
weak convergence): 

tt(/x 5 u) := inf{e > | VY G B fi(Y) < v(Y e ) + e and v(Y) < /i(F e ) + e}, 

where B is the sigma-algebra of Borel subsets of (—A, A), and Y e = U pe yS(p, e) 
is an e-neighborhood of Y. 
We claim that the map 

^4)^ j8 ([0,l) X .)6P 

is measurable; in fact, it is continuous (A small change of the coefficients 
u — (ax, 02, . . . ) in I 2 sense will yield a small change in the counting measure 
of zeroes uf in Levy - Prokhorov sense). 

Now consider the space X = P N of sequences of measures with the product 
topology. Notice that the map ft — > X defined by to \-> {v%} is measurable, 
as each coordinate is measurable; Moreover, its image lies almost surely in 
the (measurable subset) of weak converging sequences. The map C — > P 
which takes a weak converging sequence {pn) G C to its limit v G P is also 
measurable. We arrive at 

Observation 1. Any measurable set M G P induces a measurable set M = 
{co: u w G M} cCcO. 
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Consider the event: 

B = {y 6 fl : the limiting measure v w has at least one atom} cCcI 
By the last observation, B is measurable w.r.t. J 7 . 

We construct a measurable function h : B — > (— A, A) N which maps u G B 
to a list of all atoms of the limiting measure , as follows. Let h\ : B — > 
(—A, A) be the map which maps some u G B to the largest atom among 
those of (if some (finite) number of atoms share this property, return the 
left-most one). Again by observation [TJ hi is a measurable map. In a similar 
manner we construct hi, which gives the second (left-most) largest atom; and 
so forth. Our list of atoms is simply h = (hi, hi, . . . ). We notice that 

A = nj e pj{w : (^{hiu})N is a convergent sequence of numbers} =: H^Ei. 

All that remains is to prove measurability of E\. 

Indeed, the map H : X X (—A, A) — > {0, 1} which matches ({fjv}, a) to the 
indicator of the event {(vn{o>})n is a convergent sequence} is measurable; by 
composition of measurable maps 1^ = H((u%), h\uf) is a measurable function, 
as anticipated. □ 

Appendix B. Exponential Decay of Some Tail Events - 

Offord-Type Estimates 

In the course of the proof the main theorem, we used several times expo- 
nential estimates on certain probabilities: theorem[5j propositions 5.2 and 5.1 



and similar propositions for GAFs. Such estimates are sometimes referred 
to as "Offord-type large deviations estimates". We demonstrate the proof of 



proposition 5.1, as the rest are similar. We adopt the proof of Sodin [2D] , 
presented also in [5j chapter 7]. 

We first present our key-lemma, which deals with 2-dimensional Gaussian 
random variables. 

Lemma B.l. If r\ ~ A/r2 (//,£), and E is an event in the probability space 
with F(E) = p, then: 



1 pi 
-(1 + — ) \ogp + — + - \og(trace S + |/i| 2 ) 



|E(xi?log < p 

where Ai is the biggest eigenvalue ofH. 
Proof. Upper bound: by Jensen's inequality, 

-E( X £logM 2 ) <\og( E ^ 2xE ^ ) <logEy 2 -logp. 
p \ V J 

If rj — u + iv, then 

E\r]\ 2 = Eu 2 + Ev 2 = var u 4~ (Eu) 2 + var v + (Ev) 2 = trace S + 
Putting this in the previous equation, we get: 
(21) E( XE \og M) < |[log(trace S + H 2 ) - logp]. 
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Lower bound: 

E(x£logM) > -E(x B log~ \r]\) 

= -E(log~ \v\XEn{\ v \<p}) - E(log~ \v\XEn{\ v \>p}) 
The second term may be bounded below by 
(22) - E(log~ \v\XEn{\ v \> P }) > plogp 

For the first term, we begin with some general manipulations: 

[ f 1 ds 

-E(log \v\XEn{\r,\< P }) > -E(log \v\X{\v\<p}) = -E X|ij|<p / X«>|»j| — 

IP[| ?7 1 < min(p, s)] — 

Let us therefore bound from above the probability P(|^| < R). Denote by 
Ai, A 2 the eigenvalues of E, where Ai > A 2 > 0. 

r}\ < R) = j== I exp ( — —{x — /j i ) t 'E~ 1 (x — jj) J dm 2 (x) 



o 



27ry |E| J\x\<R 
< -== I exp ( — -x T Ti~ 1 x J c?m 2 (x) 

27Ta/|E| i|xi<i? V 2 J 



2tt ./| tf | <fl V 2 



12 _ K 

< / e _ 2 r rrfr = 1 — e 2A i 

o 



We have used the changes of variables y = Ux where U is the orthogonal 
matrix diagonalizing E, and later y^ = y/X~iWi for i — 1, 2. 
Continuing, we have: 

-E(log l^lxsni^Krf) > - / ds ~ I ds 

P 2 /2Ai 







1 

2 

1 /-'"' 2A ' v 1 , i, 2 , 1, 

2' 

Therefore our lower bound is 



e~*log(2Ait)cft 
> - jf log(t)dt + ^-log(2A 1 ) = ^(logp - 



E(xslog|^|) > ^( lo SP~ 2^ ~ plogp - + ^ ~ 2AT^ plogP 

Combining the two bounds we get the desired result. □ 



We now turn to the proof of proposition 5.1 Take 4>{z) = 4>t{z) a real C 
function, whose support is [— |, T+ 1] x [a', 6'] with — A <a'<a<6<6'<A, 
and which takes the value 1 on [0, T] x [a, 6). We may build such <Pt{z) that 
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will obey also the bound || A0|| L1 < 10(T + b — a). Assume a and (3 are fixed 
for now, and fix also s > 0. We are interested in dominating the probability 
of the event A T = {N T > sT}. Write p = p T = P(A T ). 
We have 

Nt<t^J A<f> T (z)teg\f(z)\dm 2 (z), 

and therefore, 

sT-p<E g ( X A T N T ) <E g [xa t ^ J A<f>(z)log\f(z)\dm 2 (z] 

= ^-J A0E 9 ( X a t \og\f(z)\)dm 2 (z) 

<i-||A0|| ia supE, ( XAt log|/(*)|) 

Before we continue, let us justify the exchange of expectation and integral. 
Recall f(z) = g(z) + T) a ^(z); so in order to use Fubini's theorem we need 
(23) 

E g \A(f)(z)-\og\g(z)+r] a:l3 (z)\\dm2(z) = / \A<j>(z)\ E g | log \g(z)+r} a>l3 {z)\ \ < oo 

Jd 

For each z G D, f(z) = g(z) + rj at p(z) is a 2 dimensional Gaussian random 
variable, with mean fi(z) = i] a ^(z), and the same covariance matrix 
as the 2 dimensional Gaussian r.v. g(z). By lemma 4.1, we see that both 
fi(z) and depend continuously on the paremeter z. So, the function 

E s | log \g(z) + 7) a ,p(z)\ is bounded above for z G support (0), which ends this 
argument. 

Notice further, that in our stationary case Xi(z), \ 2 (z), the eigenvalues of 
depend on y only, where z = x + iy. Therefore they have lower and upper 
bounds on R x [a',b']. Notice that also fi(z), being a trigonometric function, 
has such bounds. By applying lemma B.l with rj(z) = g(z) + r) a> p(z), we get: 



sup E(xa t log \g(z) + CI) < p(ci ~ c 2 logp). 

zeW.x[a',b'] 

where c\,c 2 are positive constants {c\ depending on a, (3, the horizontal lines 
a, b, and the kernel of g). Putting all this together, we get: 

5 

sT ■ p < —(T + b — a)p(ci - c 2 logp), 

7T 

which leads to the exponential bound we strived for: 

3c, C > such that p T = f g (N T > T s) < Ce~ cs , VT > 1 . 

□ 
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